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Abstract. We obtain a first order differential equation for the driving func- 
tion of the chordal Loewner differential equation in the case where the domain 
is slit by a curve which is a trajectory arc of certain quadratic differentials. In 
particular this includes the case when the curve is a path on the square, trian- 
gle or hexagonal lattice in the upper half-plane or, indeed, in any domain with 
boundary on the lattice. We also demonstrate how we use this to calculate 
the driving function numerically. Equivalent results for other variants of the 
Loewner differential equation are also obtained: Multiple slits in the chordal 
Loewner differential equation and the radial Loewner differential equation. 
The proof of our theorem uses a generalization of Schwarz-Christoffel map- 
ping to domains bounded by trajectory arcs of rotations of a given quadratic 
differential. 



Introduction 

Suppose that H — {z G C : Im(z) > 0} is the upper half-plane and 7 : [0, T") 1 — HI 
is a simple Jordan curve with 7(0) G R and 7(0, T) = {7(f) : t G (0, T)} C H. Then 
for each t G (0,T), 

H t =H\7(0,t] 

is a simply-connected domain and hence by the Riemann mapping theorem, we can 
find a conformal map f t of EI onto H t . Moreover, we can require that f t has series 
expansion 

z \z z 

Normalized in this way ft is unique and is said to be hydrodynamically normalized. 
The function C(t) is positive, continuous and strictly increasing: it is called the half- 
plane capacity of 7(0, t]. Thus we can reparameterize 7 such that C(t) = 2t for all 
t, we will call this parameterization by half-plane capacity. With this normalization 
and parameterization, the function ft satisfies the differential equation (where f' t 
denotes differentiation with respect to z and ft denotes differentiation with respect 
to t): 

(o.i) ; t( ,) = -^g) 

where = / t _1 (7(i)) 1S a continuous real-valued function. This is the chordal 
Loewner differential equation] is called the driving function of the slit 7. The 
converse is also true: given a measurable function £, the differential equation (jO.ip 
with initial condition fo(z) = z has solution f t which is a conformal map from H 
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Figure 1 . A path on the hexagonal lattice on the upper half-plane 
(left) and a plot of its driving function on the y-axis against time 
on the x-axis (right). 
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Figure 2. A path on the square lattice on the upper half-plane 
(left) and a plot of its driving function on the y-axis against time 
on the a;-axis (right). 



into itself (although /(H) is not necessarily a slit domain). Chapter 3 and 4 of [7] 
gives full details of this construction. 

Since Schramm's discovery of stochastic Loewner evolution in 1999 (see |15|). 
there has been huge interest in the chordal Loewner differential equation and its 
variants. But the relationship between the slit in HI and its resulting driving function 
is not well understood. There are a few papers that relate the behaviour of the slit 
with the behaviour of the driving function e.g. [10], [§]; also, the paper [2] calculates 
the slit arising from a few driving functions. In this paper, we will obtain a first 
order differential equation for £ (which we can then solve numerically) that allows 
us to calculate the driving function £ in the case where the curve 7 is a trajectory 
arc of a certain type of quadratic differential. We will show that this includes, for 
example, the case when 7 is a path on the square/triangle/hexagonal lattice in the 
upper half-plane or indeed, in any domain whose boundary lies on such a lattice. So 
for example, Figure Q] plots the driving function of a path on the hexagonal lattice 
in the upper half-plane and Figure [2] plots the driving function of a path on the 
square lattice in the upper half-plane. 
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We also note that we can obtain equivalent results for other variants of the 
Loewner differential equation for example, in the radial version or with multiple 
slits. We will discuss this in the paper as well. 

The proof of our formulae uses a generalization of Schwarz-Christoffel mapping 
to domains bounded by trajectory arcs of rotations of a given quadratic differential. 

We also mention that, currently, the common method used to find the driving 
function of a given slit is to use the Zipper algorithm discovered independently by 
D. E. Marshall and R. Kuhnau to approximate the function f t which can then be 
used to determine the driving function. The Zipper algorithm can be viewed as 
a discrete version of the Loewner differential equation and hence is well suited to 
studying growth processes, ft also has the advantage of being very fast. See pTj 
and [3]. 

f. Main results 

To state our main results, we have to provide some background in the theory of 
quadratic differentials. Note that not all the terms used here are standard in the 
literature. See Chapter 8 of [13] and [16] for more details. A quadratic differential 
on a domain DcC = CU{oo}is the formal expression 

Q{z)dz 2 , 

where Q(z) is a meromorphic function on D. Then for uj G D with u> ^ oo, Q(z) 
has Laurent series expansion about to, 

oo 

Q{z) = J2a k {z-uj) k 

k—n 

for some n > — oo with o„ ^ 0. Then we define the degree of lu with respect to 
Q(z)dz 2 , degg(uj), to be equal to n. 

If oo S D, then near oo, Q has Laurent series expansion given by 

oo 

Q(z) = b * z ~ k > 

k—m 

then we define the degree of oo with respect to Q(z)dz 2 , degQ(oo) to be equal to 
m — 4. The "4" in the definition ensures that the degree is conformally invariant 
in a way which we will make precise later. Then oj € D is: 

• a zero of Q(z)dz 2 if degg(o;) > 0. 

• a pole of Q(z)dz 2 if degg(o;) < 0. 

• an ordinary point of Q(z)dz 2 if degg(cij) = 0. 

A trajectory arc of Q{z)dz 2 is a curve 7 : (a, 6) 1— > D that does not meet any 
zeroes and poles of Q(z)dz 2 and satisfies 

Q{l{t))i{tf > for all t £ (a, 6). 

For 9 £ [0, 7r), a 9-trajectory arc of Q{z)dz 2 is a curve 7 : (a, b) 1— > D that satisfies 

arg[Q(7(£))7(i) 2 ] = 20 for all t € (a,b). 

Then 7 is a 0-trajectory arc of Q(z)dz 2 if and only if it is a trajectory arc of 
e~ 2lt 1 Q(z)dz 2 . Hence, a 0-trajectory arc is simply a trajectory arc and we call a 
7r/2-trajectory arc an orthogonal trajectory arc. It is clear that these definitions 
are invariant under reparameterization of 7 so we will often call the point set of 7 
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a trajectory arc or ^-trajectory arc. We call a maximal trajectory arc a trajectory 
and similarly, a maximal ^-trajectory arc is called a 6 -trajectory. For example, if we 
consider the quadratic differential ldz 2 in C, then the 0-trajectories are the straight 
lines with gradient exp(20). 

We now consider a special type of quadratic differential: Let D be a domain 
with piecewise analytic boundary. A Kiihnau quadratic differential is a quadratic 
differential, Q{z)dz 2 , on D satisfying the following two properties: 

Definition: Let D be a domain with piecewise analytic boundary. A Kiihnau qua- 
dratic differential is a quadratic differential, Q(z)dz 2 , on D satisfying the following 
two properties: 

(f ) We can write 

n 
3 = 1 

such that each Tj is an open analytic arc with F^ (1 Tj = for k ^ j 
and moreover, Q(z) extends continuously to each Tj and arg[Q(z)dz 2 ] is 
constant on each Tj i.e. each Tj is a ^-trajectory arc for some 9j G [0, n). 
(2) At z G Tj, n Tj for all j ^ k, there are either only finitely many direction 
from which trajectories approach the point z or if there are infinitely many 
directions from which trajectories approach the point z, then for each such 
direction, there is only one trajectory that approaches z at this direction. 

These quadratic differentials are studied by Kiihnau in [6] where he applies them 
to the study of certain Grotzsch-style extremal problems 

Property (i) above, also implies that D is locally connected. Thus each prime 
end of D corresponds to a unique point in dD (see [HI p. 27]). If, in addition, 
a point on dD corresponds to a unique prime end, then we make no distinction 
between the two. Let 2 be a prime end of D. Then we have 2 cases: Either z G Tj 
for some j = 1, . . . , n; or there exist exactly 2 of the (I^-), that end at the prime end 
z. In the latter case, we will denote z by Zk and assume that Tk, a ^-trajectory 
arc, and Ffc_i, a #fc_i-trajectory arc, are the only 2 arcs that end at z k . Then we 
can define the degree of z k in D with respect to Q(z)dz 2 , deg D Q(z k ), as follows: 

, 2[|0* - e k . x \/n + J fc - 1] if k + fc _i, 

deg D , Q (z k ) - j 2Jk ifd k = e k -i, 

where J k is the number of trajectories of Q(z)dz 2 inside D that end at the prime 
end Zk- If J k is infinite, then the degree is not defined. Then for prime ends z such 
that z G Tj for some j = 1, . . . , n, we define 

deg D , Q (z) = 0. 

Although the motivation for this definition currently seems unclear, we will see 
that this indeed generalizes the concept of degree to points on the boundary. In 
particular, we will show that for x G dM, if deg H q(x) G Z, then Q can be extended 
to a meromorphic function in a neighbourhood of x with 

deg H ,Q(a;) = deg Q (a;). 
We then have the following theorem on Kiihnau quadratic differentials in H: 
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Theorem 1.1. Suppose that Q(z)dz 2 is a Kiihnau quadratic differential 
Then we have 

for some constant R ^ 0, Q € C, Xj E M. for j = 1, . . . , n. 

This theorem can be viewed as a generalization of the Schwarz-Christoffel for- 
mula to domains bounded by ^-trajectory arcs of a given quadratic differential. 

We then have the following theorem on the Loewner driving function of a <fi- 
trajectory arc of a Kiihnau quadratic differential Q(z)dz 2 that starts at a point 
Co e R with deg HjQ (Co) = N e {0, 1, . . .}. 

Theorem 1.2. Suppose that Q(z)dz 2 is a Kiihnau quadratic differential on H such 
that there is a point §)6R with deg H q(£o) = N £ {0, 1, . . .}; then we have 

Q{w) = (w £ ) N (f[(w- 

where aj € C and ctj £ R. Let 7 : [0, T) t— > dM be a simple curve such that 
7(0) = £0; 7(0, T) C H and 7(0, T) is a (j) -trajectory arc of Q(z)dz 2 ((f) € [0,71-),) 
that is parameterized by half-plane capacity. Suppose that the functions ft maps 
H conformally onto H \ 7(0, t] and are hydrodynamically normalized. Then for 

te(o,r) 

(1.1) 2£(t) = -iTC-(t) -M + C+(t) - (][>i^(t)j + So, 
and 

1 • ' 5U c-(t)- m c+(t) m Mt) - m 

with initial condition £ (0) = £0 • Where the functions Aj (t) are defined by 

= ft l ( a j) f° r J = l,.--,n, 
and C + (i) > C~(t) are t/ie two preimages of £0 under f t ; 

fx ± = deg HW(0!t]! g(/ t (C ± (t))), 

and 

(n n 
fc=l / 

We can then use Theorem 11.21 to find the driving function in the case when the 
slit 7 consists of consecutive ^.-trajectory arcs of given quadratic differentials. We 
will explain how to do this in further detail later. One difficulty with using Theorem 
11.21 is that the parameterization is inherently given in terms of half-plane capacity. 
This makes it difficult to calculate the driving function £ if we do not know anything 
about the half-plane capacity of the trajectory arc (which, in general, is the case). 
The next theorem will allow us to compare the parametrization with the length of 
the slit: 
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Theorem 1.3. Suppose that Q(z)dz 2 , 7 and ft are as defined in Theorem ] l.'A Let 

0-£(*)r 

Then 7 satisfies 



(1.3) 7(*) = -24 



0(7(*)) ' 



The rest of this paper is organized as follows: In the Section 2, we will state 
some basic results from the theory of quadratic differentials and use them to prove 
Theorem ll.il Then we will use Theorem 1 1.1 1 to prove Theorems 1.2 and 1.3 in Sec- 
tion 3. In Section 4 we will discuss how to obtain the driving function numerically 
using Theorems 1.2 and 1.3. Finally in Section 5, we will discuss extensions of The- 
orem [TT^] to the case with multiple slits as well as to the radial Loewner differential 
equation. 

2. KifflNAU QUADRATIC DIFFERENTIALS AND GENERALIZED 

Schwarz-Christoffel MAPPING 



The aim of this section is to prove Theorem 11.11 We will first look at some of 
the basic results in the theory quadratic differentials that we will need. 

Transformation Law. 

Suppose that / is a conformal map from a domain D2 onto a domain D\ and 
suppose that Q\(w)dw is a quadratic differential on D\. If we define 

(2.1) Q 2 (z) = Qi(f(z))f(z) 2 

then Qi{z)dz 2 is a quadratic differential on D^. Then, it is clear that ^-trajectory 
arcs are preserved by this transformation law i.e. 

7 is a 0-trajectory arc of Qi{z)dz 2 / o 7 is a (^-trajectory arc of Q2(w)dw 2 , 

and also, for z € D-x 

de SQ 2 ( z ) = de gQ! (/(*))• 
Hence trajectories and degg are conformally invariant in the above sense. 

The following lemma tells us that the behaviour of a quadratic differential at a 
neighbourhood of a point is determined by the degree of that point. 

Lemma 2.1 (Local behaviour of quadratic differentials). Let Q(z)dz 2 be a qua- 
dratic differential on a domain D. Then for every uj G D there is a conformal 
mapping w = <fi(z) of some neighbourhood of lu such that 

dw 2 if deg Q (uj) =0, 

w n dw 2 if degg(w) = n > 1, 

Q(z)dz 2 = t w~ n dw 2 if degg(o;) = — n < —1 with n odd, 

c 2 w~ 2 dw 2 if degg(w) = —2, 

(w~ n + cw~ 1 ) 2 dw 2 if degg(w) = — n < —4 with n even. 

Here, c is the residue of a branch of Q(z) at lo. 

Proof. See Theorem 8.1 of 13 or Section 6 of [16]. □ 
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So since trajectories are conformally invariant this lemma tells us that the local 
structure of trajectories around a point u) <E D is completely determined by deggfw) 
and the converse is true as well. 

Lemma 2.2. Suppose that lj S D and degg(o;) = n. Then 

(1) For n > —1, there are exactly n + 2 trajectories of Q(z)dz 2 that end at ui 
and form equal angles with each other. 

(2) For n < —3, there are infinitely many trajectories ending at oj and more- 
over, there are \n\ — 2 directions at lo forming equal angles such that the 
trajectories approach to in these directions. 

(3) For n = —2, the behaviour depends on the value of c (as defined in Lemma 

nop . 

(a) If c is real, then the trajectories are the images of all radial lines under 
the map <f> defined in Lemma \2.1\ 

(b) // c is purely imaginary, then the trajectories are the images of all 
concentric circles under the map <f> defined in Lemma \2.1\ 

(c) // Re[c], Im[c] ^ then the trajectories are the images of logarithmic 
spirals under the map <ft defined in Lemma \2.1\ 

Proof. See Section 7 of [16]. □ 

This lemma shows that it makes sense for us to define deg^, q(x), the degree of 
a point on the boundary, in terms of the trajectories ending at a;. If D = H and Q 
extends to a meromorphic function on a neighbourhood of some x G IRU {oo} with 
deg H q(x) finite. Then by studying the trajectory structure at x, we can see that 

deg H ,Q(a;) = deg Q (x). 

This is the motivation for defining deg^j q in the way we have. The next lemma 
shows that the deg^, q is also conformally invariant: 

Lemma 2.3. Suppose that Q(z)dz 2 is a Kilhnau quadratic differential on a domain 
D and f is a conformal map of the upper half -plane H onto D. Then the quadratic 
differential Q f(w)dw 2 on H, defined by Qf(w) — Q(f(w))f'(w) 2 , is also a Kilhnau 
quadratic differential. Moreover, suppose that z € 3D is a prime end of D. Then 

degcgM = deg^Q^/^O)). 

Proof. By Caratheodory's theorem, / extends continuously to <9H and by Schwarz's 
reflection, / extends analytically across / _1 (rfe) for all k = l,...,n. Since 9- 
trajectory arcs are conformally invariant, this implies that Qf defined by (|2.1[) is 
a Kiihnau quadratic differential on H. Moreover, for all k = 1, . . . , n, f(^k) is a 
(^-trajectory arc of Qf(z)dz 2 . Also each point on dM corresponds bijectively to a 
prime end of H. Hence there is a bijective correspondence between the points of 
<9H and prime ends of D. Then 

deg D , Q (z) = deg^Q^f- 1 ^)) 

follows from the conformal invariance of trajectories. □ 

Reflection across trajectories. 

Suppose that D is a domain such that T C dD is an open interval in R. Let 
Q(z)dz 2 be a quadratic differential such that T is a trajectory arc or an orthogonal 
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trajectory arc of Q(z)dz 2 . Then let D = {z : z e D} be the reflection of D along 
T. Define 

Q~{z) = Q{z) for z G £>~. 
Then since T is a trajectory, we have 

Q(z) = Q~(z) G K for z G T. 

Thus by defining 

C Q{z) for zeD, 

(2.2) Q*(*)= < Q~{z) iorzeD-, 

{ Q(z) = Q-(z) forzer. 

it is easy to see that Q* is meromorphic in D U D~ and hence Q*{z)dz 2 is a 
quadratic differential on DUD". Thus by the transformation law (and using 
Schwarz reflection), this shows that we can extend quadratic differentials across 
trajectory arcs or orthogonal trajectory arcs. 

We will use reflection to prove the following lemma: 

Lemma 2.4. Suppose Q(z)dz 2 is a Kiihnau quadratic differential on H. Then for 
any z G cffl, deg M q(z) G TL implies that Q{z)dz 2 extends to a quadratic differential 
on a neighbourhood of z and hence 

deg HiQ (z) = deg Q (z). 

Proof. Firstly, if z G Tj for some j = 1, . . . , n. Then by definition deg M q{z) = 
and Q(z)dz 2 can be extended to a neighbourhood of z by reflection. By definition, 
every z G Tj is an ordinary point of Q{z)dz 2 and hence degg(z) = deg H g(z) = 0. 

Otherwise we write z = Zk and suppose that a #fc_i-trajectory arc, Tk-i, and a 
#/c-trajectory arc, Tk, end at Zh- Then, by definition, deg M Q(zk) € Z implies that 
Ok ~ Ok-i is a multiple of tt/2. Thus Tk-\ and T^. are trajectory arcs or orthogonal 
trajectory arcs of e~ 2t9k - 1 Q(z)dz 2 . Thus by reflection, e~ 2%dk - x Q(z)dz 2 extends to 
a neighbourhood of Zk- Hence, Q(z)dz 2 also extends to a neighbourhood of Zk- □ 

We can now prove Theorem II .11 but first, we explain briefly why we can view 
Theorem ll.f l as a generalized form of Schwarz-Christoffcl mapping: Schwarz-Christoffcl 
mapping is a method of computing the conformal map between the upper half-plane 
and a domain bounded by a polygon. See [12] for more details. If we have a confor- 
mal map / from H to some domain D such that the sides of D consist of ^-trajectory 
arcs of the quadratic differential Q(w)dw 2 . Then Q(w)dw 2 is a Kiihnau quadratic 
differential on D and hence by Lemma l2~3l Q(f(z))f'(z)dz 2 is a Kiihnau quadratic 
differential on H. Theorem II . II then implies that 

Q(f(z))f\zf = R\mz-( j )^ . 

This is precisely the Schwarz-Christoffel formula when Q(z) = 1. 

Also, we comment that the case when Q(w)dw 2 is either negative or positive on 
K (i.e. the boundary of K consists only of trajectory arcs and orthogonal trajectory 
arcs) is easy to prove: we can use reflection to extend Q(z)dz 2 to a quadratic 
differential on the Riemann sphere C. Hence Q{z) must be rational since property 
(ii) in the definition of Kiihnau quadratic differentials guarantees that Q(z) does not 
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have any essential singularities and so Q(z) is rational (since the only meromorphic 
functions on C are rational). This proves Theorem 1 1.1 1 for this case. 

Proof of Theorem \l.ll Since Q[z)dz 2 is a Kiihnau quadratic differential, we can 
find 



z\ < 



and 



(zk-i,Zk) for k = 1, ...,m, 
(z m , oo) for k = m + 1, 
(-co, zq) for k = 0, 

such that each is a ^-trajectory arc of Q(z)dz 2 for some Ok G [0, 7r). Let 

7" = {Ti, . . . , r m+ i} . 

Then take any r G 7". Since T is a (^-trajectory for some 9, T is a trajectory 
arc of e~ 2l9 Q(z)<iz 2 ; hence by reflection, we can reflect the quadratic differential 



e Q(z)dz across T to get a quadratic differential on 



{Im(z) < 0} which 



we call Q(z)dz 2 . Similarly, by rotating Q(z)dz 2 , we can reflect it across another 
T G T to get another quadratic differential Q*{z)dz 2 on H. Since Q* is obtained 
from Q by rotating twice, we have 

Q*{z) = e la Q(z) 

for some a G [0, 2tt). This shows that 

1 ' Q(z) Q*(z) 

can be extended to a meromorphic function in C\{zi, . . . , z m }. Then part (ii) of the 
definition of Kiihnau quadratic differentials implies that all the finite singularities of 
^(z) are simple poles otherwise Q(z) would have an essential singularity which, by 
the great Picard theorem, contradicts part (ii) of the definition of Kiihnau quadratic 
differentials. Thus we can write: 

*(*)=*(*)-[£; -A: 

y=i ^ 

where Q G C, and Xj G M, and h(z) is an entire function in C that does not does 
not vanish in C. This implies that 



Q(z) 



exp 



h(C)d£ 




Zk) 



Moreover, the singularity at oo of 

exp 



h(C)dC 



cannot be essential by part (ii) of the definition of Kiihnau quadratic differentials 
(otherwise we would get a contradiction with the great Picard theorem as above). 
This implies that 



exp 



h(()d( 
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is constant (since it has no zeroes or poles). Hence 

□ 

If £j € C \ R, then by definition, we must have 

*j = degg(Ci). 
Moreover, if Q £ R and deg W g(Cj) < oo we also have 

v k = deg mQ {z k ). 

We will not prove this fact here but in the following corollary we will consider a 
special case. The general proof follows readily from it. We will prove the following 
corollary which is simply an application of Theorem 11.11 to domains slit by (j>- 
trajectory arcs: 

Corollary 2.5. Suppose that Q{z)dz 2 is a Kiihnau quadratic differential on H such 
that there is a point (oeR with deg H g(^o) = N 6 {0, 1, . . .}; then we can write 

(2.3) Q(w)=R(w-So) N \f[(w 

where dj € C, ctj £ R, and also R is some non-zero constant. Let 7 : [0, T] 1— » H be a 
simple curve such that 7(0) = £0 and 7(0, T) is a <p -trajectory arc of Q(w)dw 2 in H 
(4> S [0, 7r) ) and £ = 7(T) G H is an ordinary point of Q(w)dw 2 (i.e. degg(C) = 0). 
Suppose that f maps H conformally onto H\7(0,T]. Then f satisfies 

n 

(2.4) Q(f(z))f'(z) 2 = R'{z - 0\z ~ c~r~ (z - c+Y + JJ(z - , 

i=i 

where R' is some constant; c~ , c + are the two preimages of £0 under f satisfying 
c~ < c + ; Aj is £/ie preimage of aj under f ; and £ is t/ie preimage of £; and 

= d eg H \ 7 (o, s ],Q(/( c± ))- 

Proof. Theorem II .11 and Lemma [2~4l imply that Q(w) can be written as (|2.3p . Then 
Lemma T2.3I implies that Q(z) = Q{f{z))f'{z) 2 is a Kiihnau quadratic differential. 
So by Theorem II. 1[ we only need to look at the singularities of Q(z)dz 2 . 

Now, by Schwarz reflection, / extends to a conformal map on C \ {£, c _ ,c + }. 
Thus, by the conformal invariance of trajectories, this implies that 

n 

Q f (z) = R'(z - M (z ~ c-y~{z - c+Y + ]J(z - A,-) ". 

3=1 

Then by Lemma 12.21 there are exactly two ^-trajectory arcs of Q(z)dz 2 ending 
at C = l{T) of which 7(0, T) is one of them. So by the conformal invariance of 
trajectories, there is one ^-trajectory arcs of Q(z)dz 2 ending at £ that is contained 
in H. Hence, by definition, degg- H (£) = 2. Using Lemma [2.41 this implies that 
degg(£) = 2 i.e. M = 2. Thus we only need to determine and [i + . 
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Note that since £o has degree N with respect to Q(z)dz 2 , we can determine, 
using Lemma l2~2"1 that the angle between 7(0, T) and /((£, c - )) at £0 is 

deg Hi g(/(c-)) + 2 

TV = T ~t , 

y N + 2 ' 

and similarly, the angle between 7(0, T) and /((£, c - )) at £0 is 

de gifiQ (/( C +)) + 2 

7rt/» = 7T . 

r N + 2 

Hence, by Schwarz reflection, the function 

F(z) = (f(z)-^~ 

extends to a conformal mapping on a neighbourhood of c~. Thus in a neighbour- 
hood of z = cT , we can write 

(2.5) f{z) = Z + {z-c-)^h{zf- , 

where /i is analytic in a neighbourhood of z — c~ with h(c~) =/= 0. Now 
Q' f (z) Q'(f(z))f'(z) 2 , /"(*) 



Q/(*) Q(/O0) /'(*) ' 

The residue at z = c~ of the left-hand side of the equation is fi~ , and we can use 
to determine the residue at z = cT of the right-hand side. Thus we get 



V- = dc EH,Q(f( c ))• 

We apply the same method to c + to get fi + . □ 

3. Domains slit by ^-trajectory arcs 

Let Q(w)dw 2 be a Kuhnau quadratic differential on H with deg H q(£o) = N £ 
{0, 1, . • •} for some £0 G Then by Theorem II. II and Lemma[ 



Q(w) = (w- Co) 



TV 




where a>j £ R and dj 6 C. Now suppose that 7 : [0, T) 1— > H is a simple curve 
such that 7(0) = £0 and 7(0, T) is a </3-trajectory of Q{z)dz 2 in H (0 = [0,7r)) 
that is parameterized by half-plane capacity. As mentioned in the introduction, 
there exists conformal maps ft ■ H 1— * Ht = H \ 7(0, t] satisfying the hydrodynamic 
normalization. Then by restricting Q(w)dw 2 to a quadratic differential on H t we 
can induce via f t and (|2.1[) , a quadratic differential on H: 

(3.1) Qt(z)dz 2 = Q{f t {z))f' t {z) 2 dz 2 . 

We now use Corollary 12.51 and (|3.ip to prove Theorem 11.21 



Proof of Theorem Note that by Schwarz reflection, each f t can be extended 
to a conformal map on C \ {C~(t), C + (t), £(i)}. Then since ft{z) satisfies the 
hydrodynamic normalization, this implies that 



m 2 = 1+0 (^) 



00. 
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So by (JSU), 

(3.2) jftMj.i + o^)..,^. 

If we let C = then we get 

Since ft is analytic in a neighbourhood of infinity (|3.3[) is a Taylor series expansion 
and hence we can look at the Taylor series coefficients, in particular: 



C(0,e) C 2 * y C(0 , e) C 2 Q(/t(l/0) 

for small enough e > where C(0, e) is the anticlockwise contour about the circle 
with centre at zero and radius e > 0. Then by Theorem 11.11 and Corollary 12.51 we 
can write 

Q(w) = R(w - £ ) N (Ww- aj )A , 

and 

n 

Q t (z) = R'(z - if{z - C-(t)T (z C+(t)Y + l[(z Aj(t)) aj . 

j=i 

Hence by the residue theorem (since /t(l/C) = 1/C + ' ' ' as C ~~ ¥ 0), this implies 
that 



fc=i l\ L \fc=i / \;=i 



= o. 



This implies (II. To get (| 1 . 2[) . note that / t satisfies the chordal Loewner differen- 
tial equation (jO.ip and hence if we let Gt — ft' 1 ° f s f° r some s G (0,T) fixed and 
t > s, then the chain rule implies that Gt satisfies the differential equation 

dG tf . 2 



dt v ' G t (z)-£(t) 

Then for some s sufficiently close to t, we can write each Aj(t) = Gt(wj) Wj G C 
for all j = 1, ... ,n. Thus 

Aj{t) = a^twy 

Similarly we get 

d±{t) = c±$TWy 

Hence we get (| 1 . 2|) from differentiating (jl.ip . □ 
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An extension: 

We can extend Theorem ll.2l to the case when 7 is made up of different ^-trajectory 
arcs of some quadratic differential Q(z)dz 2 : Let 7 : (0, T] 1— > H be a curve with 
7(0) € R such that there is a partition 

{0 = t < h < ■ ■ ■ < t r = T} 

such that 7(tfc_i,ifc) is a ^-trajectory arc of Q(z)dz 2 and j(tk) is an ordinary 
point of Q(z)dz 2 for k = 1, . . . , r. Then we can find the driving function £(t) of 7 
by applying Theorem 11.21 to the #i-trajectory arc 7(0, t\) to get a driving function 
£i(i), and applying Theorem 11.21 inductively to each fZ (^(tk, ife+i)) (which is a 
0fc+i-trajectory arc of the quadratic differential Qt k (z)dz 2 = Q(ft k (z))f^ (z) 2 dz 2 ) 
to get £ k (t). Then 

£(*) = &(*) forte [* fc _i,tjfe). 
We also have the following corollary: 

Corollary 3.1. Suppose that Q(w)dw 2 and 7 are as defined in Theorem 
Then the driving function £ and Aj , C~ , C + as defined in Theorem are in 
C°°(0,T). Moreover, we can write any derivative of £, C~ , C + , Aj explicitly in 
terms of £, C™, C + ,Aj and the exponents aj . 

Proof. Recall that, in the proof of Theorem 11.21 we had the formulae 



This implies that each term in (|1.2p is diffcrentiable so we can write the second 
derivative of £ in terms of £(t),Aj(i), C^{t). This in turn implies that we can write 
the third derivative of £ in terms of £ (i) , Aj (t) , (t) and the exponents. Continuing 
inductively, we have showed that every derivative of £ exists and can be expressed 
in terms of Aj(t), C ± (i) and the exponents. Note that each derivative of £ is 
finite for t G (0, T) since 

\A j (t)-mi\c ± (t)-m\>0' 

Then £ is smooth implies that Aj(t), C ± {f) are also smooth. □ 
Theorem 1 1 . 31 then follows from Corollary 12.51 and Theorem 11.21 



Proof of Theorem \1.3[ 

First note that, by the definition of ^-trajectory arcs, 7 always exists and is never 
0. Also by Corollary |2~51 <& t (£(*)) ^ 0, 00; thus the right hand side of fO]) always 
exists since, by definition, 7 avoids poles and zeroes of Q(w)dw 2 . 
Recall that /t(£(t)) — j(t), this implies that 

at) = Mm +nm)m 

Then combining the Loewner differential equation (|0.ip with (|3.ip we have 



/t(*) = - 



Qt(z) **(*) 



z-m\Q(ft{z)) y Q(f t (z)) 



Mam = -2a 



. Q(7(<)) 

Note that $t(£(t)) ^ 0,oo since, by Corollary |2.51 Qt(.z) has a double zero at 
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Thus we have 



7(<) = 2 v + v qwo) ?(<) - 2 v w 

since Theorem 11.21 implies that £ is finite for all t G (0,T) and Corollary 12 . 51 implies 
that Qt(f(*)) =0. □ 

4. Applying Theorem 11.21 



In practice, understanding £(t) via is not possible: it is difficult to calculate 
the positions of the zeroes and poles of Qt because the information we have on them 
is all relative to £(t) (which we are trying to find). On the other hand, (|1.2p is more 
useful in applications. In this section, we will demonstrate how we can use (ll.2[) 
to calculate numerically the driving function of a given slit that consists of 9 k - 
trajectory arcs of a given quadratic differential. The method is basically a modified 
version of Euler's method. 

Firstly, for any smooth function h on (0,T), Taylor's theorem implies that for all 
M =1,2,..., 
(4.1) 



h 1 1 + \- ) - I Mt) + > ; —— - i < — — 7 sup 



K 




se(M+i) 



d M h, . 



for t,t + 1/K G (0, T). We will apply (|4~Tj) to the functions £,A k and C ± (as 
defined in Theorem 1 1.2[) noting that, by Corollary 13. 1( they are smooth and all of 
their derivatives can be expressed in terms of £(t),Ak(t) and C ± (i). Thus if we 
know £(s), A k (s) and C (a) we can use (|4.1[) to obtain an approximate formula for 
£(s + A -1 ), A fe (s + if -1 ) and C ± (s + A" 1 ) (choosing A to be small and/or M to 
be large so that the right-hand-side of (14. lj) is small); then we can apply (|4.1| to 
^(s + if- 1 ), A fc (s + A" 1 ) and C^s + A" 1 ) to find £(s + 2A" 1 ), A fc (s + 2A" 1 ) and 
C ± (s + 2A _1 ). Continuing like this, we obtain an approximation of £ at the points 
{s + nK- 1 }. 

So clearly what we need to do now is find the starting values £(s),Ak(s) and 
C ± (s) so we can apply the above method. But because £ is not differentiable at 
0, we cannot use the formula (|4.ip with t — 0. The way around this is to note 
that if deg H q(£o) =JVe {0, 1, 2, . . .} then since we know deg H q ( (C + (t j), we can 
calculate the angle that the trajectory makes with the line [£o, oo) (as in the proof 
of Corollary |2.5p . Then we find that the angle is 7r0 where: 



TCIp — 7T 



2deg HjQt (C+(i)) + 2 



N + 2 

So if we choose s small enough, we have 

where Ff^° is the conformal map that maps H conformally onto Hf that is hydro- 
dynamically normalized where Hf is the upper half-plane slit by the straight line 
starting at £o making an angle trip with [£o, oo), with half-plane capacity 2s. Then 
we also have 
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and also, C~(s),C + (s) are approximately the two preimages of £o under Ff + '^°. 
Then we can use (jl.ip to calculate £(s) approximately. We can then plug this 
information into (|4.ip as described above. 

Note that Ff' x can be found using the fact that 

(4.2) Flf{z) = {z-(l-2p)^t- py/tf (z - (1 - 2p) Vt + (1 - p)y/t) ^ 

for some A. Then we reparameterize this formula to remove the A and translate the 
point to £o- Unfortunately inverting this function cannot be done explicitly but 
it can be done numerically very efficiently using Newton's method. Alternatively, 
by selecting a small s, we can assume that 

A k (s) « a k 

for all k. Then we note that the 2 preimages of £o under Ff' x can be determined 
explicitly (see jTTj). This obviates the need to numerically invert Ff' x . 

Another difficulty is that, in general, given a slit, we cannot parameterize it by 
half-plane capacity so it would be difficult, for example, to know at which t one 
should stop. Most formulae for calculating half-plane capacity of some compact set 
K rely on knowing the conformal map fx of H onto M\K (normalized hydrodynam- 
ically). One possibility would be to use the probabilistic definitions of half-plane 
capacity given in [7]. We will use the fact that Theorem 1 1.31 and Corollary [3J] imply 
that we can give all derivatives of j(t) in terms of Ak(t), C~(t), C + (t) and the 
exponents ^,p + ,a k so if we know these, we can also use (|4.1I) to approximate 
7. This in turn allows us to calculate the length of the slit 7. Thus if we know 
beforehand length of our slit, we can calculate at what value of t we stop. 

We now have everything we need in order to use (|4.1|) to calculate the driving 
function numerically of any slit that is made up of Q k trajectory arcs of a qua- 
dratic differential Q{w)dw 2 . We will demonstrate how this is done in the following 
example: 

An example. Suppose that 7 : (0, T) — > EI is a piecewise linear arc parameterized 
by half-plane capacity that satisfies: 

• 7(0) = 0. 

• From t = to t = t±, 7 is the straight line arc from to i; call this Fi. 

• From t = t\ to t = 7 is the straight line arc from i to 2 + i; call this T2. 

• From t = ti to t = t% = T, 7 is the straight line arc from 2 + i to 2 + 2i; 
call this 

First note that 7 is made up of alternating (jr/2)- and 0-trajectory arcs of the qua- 
dratic differential ldw 2 in EI and hence we can use Theorem ll.2l (or more specifically 
the extension of Theorem 11.21 detailed in Section 3) to calculate £. As mentioned 
previously, there is no easy way to know beforehand what t\ , . . . , t3 are. For sim- 
plicity, we will only use M = 1 in (|4.1jl i.e. 



and fix a large K . Obviously Ti forms a right angle with real line; so we can use 
(14.21) to determine the function 



f tl = F?/ 2 '°(z) = v^^r. 
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-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



Figure 3. The example path in the upper half-plane (left) and a 
plot of its driving function on the y-axis against time on the x-axis 
(right). 

It is easy to see that in this case, t\ = 1/4 and £ is constantly for t £ (0,ii]. 
This induces the quadratic differential using (|3.1[) : 

z 2 dz 2 

Q tl (z)dz 2 = — . 

(z + l)(z-l) 

Hence, we let Ai(ti) = —1, A 2 (ti) — 1. Also / 4 ~ (72) is a 0-trajectory arc of 
Qi(z)dz 2 starting from £(£i) = onl (by the conformal invariance of trajectories). 
Now note that / t 7 1 (72) makes an angle of 7r/4 with the positive real axis, and hence 

,1/4,0 



since A is large. We can then use Newton's method to find the preimages under the 
above approximation of f tl +K- x °f the points A\(ti), ^2(^1) and the 2 preimages of 
zero to get the points A^ti + K A 2 (t 1 + A" 1 ), C~(h + A'" 1 ), C + {h + A" 1 ) and 
hence, using we can find ^(ii+A -1 ). Then inserting this into (|4.ip . as detailed 
above we can also find + nK^ 1 ) and A\(t\ + nA _1 ), A 2 (ti + nK^ 1 ), C~(ti + 
nK^ 1 ), C + (ti + nA _1 ); also, by Theorem ll.3[ we can find |7(ti + nA _1 )| if we let 



1 



3=1 



-nK- 1 ^ > (length of T 2 ) = 2 



then t 2 (K) as t 2 for A large. So we just assume that t 2 — t 2 (K). Let A^(t 2 ) = 
C-{t 2 ) and A 4 (t 2 ) = C+(t 2 ). Hence by (I3~TD . 

n , w 2 _ (z-^2)) 2 (^-A 3 ft2))^ 2 

Vt21 j (^ - iMtaJX* - A a (ta))(a - A 4 (t 2 )) - 
Then, by the conformal invariance of trajectories, Q?a) is a 7r/2-trajectory of 
Qt 2 (z)dz 2 and also, /T (r 3 ), makes an angle 37r/4 with (£(t 2 ),oo) and so 

Then, as before, we can use Newton's method to find the preimages under the above 
approximation of f^+K- 1 of the points Ai(t 2 ), . . . ,Ai(t 2 ) and the 2 preimages of 
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£(t 2 ) to get the points A^+R' 1 ), . . . , Aa^+R' 1 ), C'fa+R- 1 ), C+fa+R- 1 ) 
and hence use (jl.ip to get £(£2 + We insert these into the formula iteratively 

to get Cfo+nK- 1 ) and A^+nR- 1 ), . . . , Aifa+nR- 1 ), C~ (t 2 +nR- 1 ),C+(t 2 + 
nR^ 1 ) until t 2 + nR~ 1 « T. Thus the end result is that we found the driving func- 
tion of the first 3 steps of the slit given in Figure [3] Of course, our calculation of £ 
will be more accurate by taking larger R. 



For example, we can use the above method to calculate the driving function of 
any path on the square/trianglc/hcxagonal lattice on H starting from some point in 
R. In fact we can calculate the driving function of a path on the square/triangle/hexagonal 
lattice in any polygon D by mapping the half-plane conformally onto D and pulling 
back the quadratic differential ldw 2 on D to Q(z)dz 2 on HI using the transforma- 
tion law. Also note that, in general, any curve 7 can be approximated by a curve 
75 which lies on the square lattice 51 2 . Then it can be shown that 

£5 — ► £ uniformly as 5 \ 0, 

where £5 is the driving function of 75 and £ is the driving function of 7 hence, we 
can use the above method to calculate £5 then take the limit as S \ to obtain £. 

Another point to note is that using the above method, we do not need to know 
before hand what the trajectory arc of the given quadratic differential looks like; 
so for arbitrary Kiihnau quadratic differentials, we can use this method to plot the 
trajectories starting at the boundary. 

We end this section by looking at what happens when the slit approaches the 
boundary: 

Proposition 4.1. Suppose that 7 : [0, T) i— > MI is a simple curve such that 7(0) € K 
and 7(0, T) is a Q -trajectory arc of some quadratic differential Q(z)dz 2 . Then let £ 
be the driving function of 7. If 

lim 7 (<)eIRU7(0,r), 

IT J 



i.e. 7 makes a loop at time T. Then 



d n f 
d n V 1 



00 as t y T 



for all n = 0,1, . 



Proof. 

For t £ (0,T), we define 

T(t)=Ws):s€(t,T)}. 

Then T t is a (^-trajectory arc in H t = H\7(0, t] of Q(w)dw 2 and it is also a crosscut 
in H t (see [H]). Then by the conformal invariance of 6>-trajectories, / t _1 (r t ) C H 
is a 0-trajectory arc of Q t (z)dz 2 . Moreover, / t _1 (r t ) is a crosscut of H with one 
end point at and the other end point in R such that either C + (t) or C~(t) 
is contained in the closure of the bounded component of H \ / t _1 (r t ). Without 
loss of generality, assume it is C + {t). Then since diam(/ t _1 (r t )) -> as t / T, 
we must have £(t) = C+(T) and hence by -> 00 as t / T. Similarly, 

we differentiate (|1.2[) as mentioned in Corollary 13.11 to obtain the result for higher 
order derivatives. □ 
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This means that as 7 gets closer and closer to making a loop, the approximation 
by (|4.1j) stops working no matter what M we choose. This phenomenon can be 
observed in Figure 02 as we turn the last corner in 7, we can see that £ decreases 
faster even though the slit is not yet that close to the boundary. 



5. Generalizing Theorem 11.21 

5.1. Multiple slits. Suppose that jk ■ [0, T) — ► H for A; = 1, . . . , N are disjoint 
simple curves such that 7fe(0) £ K and 7/5(0, T) c H. By the Riemann mapping 
theorem, there exists unique f t that map H conformally onto H t = H\UfeLi 7fe(0> t] 
that satisfies the hydrodynamic normalization. We can reparameterize such that 

N 
k=l 

has half-plane capacity 2t. Then ft satisfies 



N 



(5.1) 
where 

k=l 

and £k(t) = ft~ 1 (jk(t))- See [5] for more details. 

Theorem 5.1. Suppose that Q(w)dw 2 is a Kilhnau quadratic differential on 
such that the points £fc(0) £ R satisfy 

deg H)Q (£fc(0))=/? fc e {0,1,2,...} 

for all k. Then we can write 

\fe=l / \j=l 

with dj S C and ctj £ M. Then suppose that 7^ : [0, T) — > H /or fc = 1, . . . , iV are 
disjoint simple curves such that 7fe(0) S K and 7fe(0, T) C H and ore parameterized 
as above. Then 

N I N \ / n \ 

(5.2) 2^> = - - £a^(*) + E , 

fe=i Vfe=i / \j=i / 



and 
(5.3) 



ft (t) - ft (*) J 2 ^ <?- (t) - 6 (f ) C+ (t) - (t) J 2 Aj (t) - ft (t) 



/or all I £ {1,...,N}. Where C k (t) and C k (t) are the two preimages o/£fc(0) 
under f t satisfying C^(t) < Cf{t); 

^=deg HuQ (f t (Ct(t))); 
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Aj(t) = f t 1 (a j ); and 

Proof. By Theorem 11.11 and Lemma \2A\ we can write 

Q( W ) = (j[(w - Ckm^ (f[(w- aj r^ 

Then either by modifying the proof of Corollary 12.51 or iterating N slit functions 
and applying Corollary 12.51 N times, it is not too difficult to see that if we define 
Qt(z) by (J3T]), then 
(5.4) 

Vfc=i / y=i 

Then the proof of ()5.2p is exactly the same as the proof of (|1.1[) in Theorem 11.21 
The proof of (|5.3j) is more complicated. First let 

— am 

Then (|5.ip becomes 

/,(*) = fl(z)Pt(z). 

Now take the logarithmic derivative of Qt(z) with respect to z and t separately 
using the definition of Qt(z) given by (13. ip to get 

Q't(z) _ Q'(ft(zM(z) fj'jz) 

Q t (z) Q(f t (zj) + f'{ Z y 

and 

Qt(z) _ Q'(Mz))f t (z) +2 K(z) 



Q t (z) Q{f t {z)) fl(z) 

Q'(Mz))fl(z)P t (z) n>(z)P t (z) + Mz)P{(z) 
Q(ft{z)) fl(z) 

where we substitute (|5.ip in for f t to get from the first to the second line. Thus we 
get 

/r Qt(z) Q't( z ) u / ^ , od'c ^ 

(5 ' 5) Q^)^Q^) Pt{z) + 2Pt{z) - 

So then by (|5.4p . we note that 

0*^) ~ \h z - + z-C^(t) + z- C+(t) ) + \f^'z- Mt) j ' 

and 

Qt(z) 
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Figure 4. A plot of the two driving functions £1 (top) and £2 
(bottom) on the y-axis against time on the x-axis. 



Thus substituting this into (|5.5[) and comparing the coefficient of 

1 

(i.e. the residue at z = £i(t) of both sides of (|5.5[) L we find that this is exactly 
(IOD . □ 

Similarly, we can prove a version of Theorem II .31 and Corollary 13.11 for multiple 
slits. This means that we can use the method detailed in Section 4 with (|5.3[) 
to calculate the driving function for multiple ^.-trajectory arc slits. For example 
Figure [4] plots the graph of the two driving functions £1 and £2 m the case when 
71 and 72 are 2 vertical slits starting from -1 and 1 (i.e. orthogonal trajectories of 
ldz 2 ) and growing at the same speed. Compare this with Figure 7 in [2]. 



5.2. Radial Loewner evolution. The chordal Loewner differential equation was 
introduced because the upper half-plane was an easier domain to work with for 
many applications but the original setting of the Loewner differential equation is in 
the unit disc D = {z € C : \z\ < 1}: Suppose that 7 : [0, T) 1— > D is a simple curve 
such that 7(0) e T = {z : \z\ = 1} and 7(0, T) C B \ {0}. Then D t =D\ 7 (0, t] 
is simply-connected and € D t for all t £ (0,T). Hence the Riemann mapping 
theorem implies that there is unique conformal map ft mapping D conformally onto 
D t such that /t(0) = and / t '(0) > 0. Then Schwarz's lemma and the Caratheodory 
kernel theorem implies that / t '(0). is strictly decreasing and continuous so we can 
reparameterize such that //(0) = e~*. / t '(0) is sometimes called the conformal 
radius of D t ; hence in this case we are parameterizing by conformal radius. Then 
the functions ft satisfy the radial Loewner differential equation: 



See [9] for more details. 
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Theorem 5.2. Suppose that Q(w)dw 2 is a Kiihnau quadratic differential on D such 
that deg Q (0) = K e 1 and e^° eT = {|z| = l} satisfies 

deg D , Q (e«°) =Ne {0,1,2,...} 

then we have 

Q( w ) = w K (w- e ie °) w f f[(w - a,-)"' 

V=i 

where dj <E C and ctj € K. T/ien £/ 7 : [0,T) 1— * D is a simple curve such that 
7(0) = and 7(0, T)cD \ {0} is a <fi -trajectory arc of Q(w)dw 2 in D that does 
not meet and is parameterized as above. Then we have 

n 

(5.6) e 2 ^ (t) = e~ 2t n C- (t)-^ C + (t)-" + YlAjit)-^, 

and 
(5.7) 

• /N I ( C-(t) + e*W ,C7+(t)+e*W ^ A,(t)+e*«*> 



C-(t)-e«(*) r C+(t)-e«W ^ ^(t)-e««W 
where, as usual, the functions Aj (t) are defined by 

A i(t) = /i _1 ( a j) for j = n, 

and 

H ± =deg DuQ (f t (C ± (t))), 
C + (t) > C~ (t) are the two preimages of under f t ; and also, 

n 

Il =e iN t°l[a?. 

i=i 

Proof. The formula for Q(w) can be obtained from Theorem 1 1.1 1 bv the transforma- 
tion law. We then define Q t by (|3.1|) . Since the point is fixed by f t , this implies 
that the degg t (0) = degg(O) = K. Thus we can apply Corollary [23] (again, using 
the transformation law) to get 

Q t (z) = z K (z - e*Wf(z - C-{t)Y~ (z - C+{t)Y + \J[(z - A 3 (t))^ 

Then since by definition, 

Q t (z) = Q(f t (zM(z) 2 . 

This immediately implies (|5.6p by substituting z = 0. Then we get (|5 . T[) in the 
same way as we get (| 1 . 2|) from (jl.ll) in the proof of Theorem 11.21 □ 

As in the case of multiple slits, a version of Theorem 11.31 and Corollary 13 . 1 1 holds 
for this case. 
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5.3. Other versions of the Loewner differential equation. There are several 
other versions of the Loewner differential equation for simply-connected domains in 
the literature; the methods in this paper should work in those cases as well and the 
proofs should be similar to the proofs of Theorem [L2] etc. Also, [4], [5] generalizes 
the Loewner differential equation to multiply-connected domains and again, some 
of the methods should work in these cases possibly using methods in [Ij to extend 
Theorem 11.11 to multiply-connected domains. Finally, even if we consider general 
2-dimensional growth processes given by the Loewner-Kufarev differential equation 
(see Chapter 6 of [13]), some of the methods in this paper should still be applicable. 

Acknowledgement. The author would like to express his gratitude to his super- 
visor Dr. T. K. Carne for his constant guidance and helpful discussion when writing 
this paper. 
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